home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / dlasq1.z / dlasq1
Encoding:
Text File  |  2002-10-03  |  3.7 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAASSSSQQQQ1111((((3333SSSS))))                                                          DDDDLLLLAAAASSSSQQQQ1111((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLASQ1 - compute the singular values of a real N-by-N bidiagonal matrix
  10.      with diagonal D and off-diagonal E
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DLASQ1( N, D, E, WORK, INFO )
  14.  
  15.          INTEGER        INFO, N
  16.  
  17.          DOUBLE         PRECISION D( * ), E( * ), WORK( * )
  18.  
  19. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  20.      These routines are part of the SCSL Scientific Library and can be loaded
  21.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  22.      directs the linker to use the multi-processor version of the library.
  23.  
  24.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  25.      4 bytes (32 bits). Another version of SCSL is available in which integers
  26.      are 8 bytes (64 bits).  This version allows the user access to larger
  27.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  28.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  29.      only one of the two versions; 4-byte integer and 8-byte integer library
  30.      calls cannot be mixed.
  31.  
  32. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  33.      DLASQ1 computes the singular values of a real N-by-N bidiagonal matrix
  34.      with diagonal D and off-diagonal E. The singular values are computed to
  35.      high relative accuracy, in the absence of denormalization, underflow and
  36.      overflow. The algorithm was first presented in
  37.  
  38.      "Accurate singular values and differential qd algorithms" by K. V.
  39.      Fernando and B. N. Parlett, Numer. Math., Vol-67, No. 2, pp. 191-230,
  40.      1994,
  41.  
  42.      and the present implementation is described in "An implementation of the
  43.      dqds Algorithm (Positive Case)", LAPACK Working Note.
  44.  
  45.  
  46. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  47.      N     (input) INTEGER
  48.            The number of rows and columns in the matrix. N >= 0.
  49.  
  50.      D     (input/output) DOUBLE PRECISION array, dimension (N)
  51.            On entry, D contains the diagonal elements of the bidiagonal matrix
  52.            whose SVD is desired. On normal exit, D contains the singular
  53.            values in decreasing order.
  54.  
  55.      E     (input/output) DOUBLE PRECISION array, dimension (N)
  56.            On entry, elements E(1:N-1) contain the off-diagonal elements of
  57.            the bidiagonal matrix whose SVD is desired.  On exit, E is
  58.            overwritten.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDLLLLAAAASSSSQQQQ1111((((3333SSSS))))                                                          DDDDLLLLAAAASSSSQQQQ1111((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      WORK  (workspace) DOUBLE PRECISION array, dimension (4*N)
  75.  
  76.      INFO  (output) INTEGER
  77.            = 0: successful exit
  78.            < 0: if INFO = -i, the i-th argument had an illegal value
  79.            > 0: the algorithm failed = 1, a split was marked by a positive
  80.            value in E = 2, current block of Z not diagonalized after 30*N
  81.            iterations (in inner while loop) = 3, termination criterion of
  82.            outer while loop not met (program created more than N unreduced
  83.            blocks)
  84.  
  85. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  86.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  87.  
  88.      This man page is available only online.
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.